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Abstract The complex shear modulus of solid '^He exhibits an anomaly in the 
same temperature region where torsion oscillators show a change in period. We 

Qj • propose that the observed stiffening of the shear modulus with decreasing temper- 

ature can be well described by the response of glassy components inside of solid 
^He. Since glass is an anelastic material, we utilize the viscoelastic approach to 
describe its dynamics. The viscoelastic component possesses an increasing relax- 
ation as temperature decreases. The response functions thus derived are identical 

^, • to those obtained for a glassy, time-delayed restoring back-action. By generaliz- 

ing the viscoelastic equations for stress and strain to a multiphase system of con- 
stituents, composed of patches with different damping and relaxation properties, 
we predict that the maximum change of the magnitude of the shear modulus and 
r^ , the maximum height of the dissipation peak are independent of an applied external 

PJ ' frequency. The same response expressions allow us to calculate the temperature 

O , dependence of the shear modulus' amplitude and dissipation. Finally, we demon- 

, ^, ■ strate that a Vogel-Fulcher-Tammann (VFT) relaxation time is in agreement with 

available experimental data. 
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1 Introduction 

The low-temperature anomaly of solid helium reported for torsional oscillator 
(TO) by Kim and Chan ''^ is often regarded as evidence for supersolidity. In addi- 
tion, the many experiments that followed confirmed that defects in solid "'He are 
important to produce an anomaly. Direct experimental evidence for a true phase 
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j^ ' -Center for Integrated Nanotechnologies, Los Alamos National Laboratory, Los Alamos, New 

Mexico 87545, USA 



transition into a supersolid state, on the other hand, remains inconclusive. To date 
no definitive sign of Bose -Einstein condensation (BEC) has been seen in measure- 
ments of the mass flow^*^, the meking curve ^, and the lattice structure. ^-^ 

Controversies exist in interpreting the anomalies, since defects can display 
their own dynamics and contribute to observables in the parameter range where 
supersolidity is expected. A detailed analysis on the behavior of defects is thus 
needed to identify the existence of a supersolid state. We propose a theoretical 
framework based on a possible glass component^'^*' to capture the dynamics of 
defects in solid ^He'^'^^. This glassy component is suggested to be the cause for 
the TO and thermodynamic anomalies reported so far. Further, it is consistent 
with reported signatures of long equilibration times, hysteresis, and a strong de- 
pendence on growth history. We demonstrated in previous worb^iifliiiidii^ that the 
freezing out of defect dynamics can account for the anomalies in thermodynamic 
and mechanical experiments. The mechanical experiments like torsion oscillator 
and shear modulus were treated by introducing a back-action term with time delay. 

Based on the same glass hypothesis, we use the viscoelastic approach in this 
paper to describe the latest shear modulus experimental results. The usage of vis- 
coelastic approach is justified since glass is a viscoelastic material. In fact, our 
approach here is similar to the model proposed by Yoo and Dorsey '^ for the TO 
experiments. We propose the presence of a distribution of viscous components 
embedded in an otherwise elastic solid. The model we introduce is also known 
as the generalized Maxwell model. It leads to the same shear modulus expres- 
sion as the one previously derived within the back-action formulation^^. This is 
sensible since the back-action describes how the glass, which is an anelastic ma- 
terial, responds to the external shear stress. In addition, we observe that (OT is 
the scaling parameter of the response functions and we discuss the consequences. 
Finally, we compare our calculations with available measurements of the shear 
modulus.iS.i2^ 



2 Viscoelastic Model 

The anomalous low-temperature stiffening of the shear modulus can be obtained 
from a viscoelastic approach. The equivalent circuit model is sketched in Fig. [T] 
We assume that the glassy components of the solid give rise to a viscous con- 
tribution on top of the dominant elastic behavior of crystalline helium. Here we 
describe the coupling between the glassy subcomponent and solid helium by a 
generalized Maxwell model, shown in Fig. |2l For pedagogical reasons, we start 
by considering a single Debye relaxor as plotted in Fig.[TJa). The Debye relaxor 
is composed of a serial connection of a rigid solid (RS) part, characterized by an 
elastic shear modulus jUrs, and a Newtonian liquid (NL) part, characterized by a 
viscosity T] . The RS part describes the ideal elastic solid helium, while NL repre- 
sents the glassy component, which gives rise to viscous damping. The two parts 
are connected in series, so that both share the same magnitude of stress, while the 
net strain is additive. The strain rate equation for both constituents is 

£ = a/iJ.Rs + a/ri, (1) 

where £ is the net strain of the Debye relaxor and a is the magnitude of stress 
shared by the components RS and NL. In order to obtain the above equation, we 
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Fig. 1 Sketch of a Debye relaxor. (a) A single Debye relaxor element composed of an elastic 
spring (/Xrs) in series with a damped dash-pot (tj). (b) A single Debye relaxor connected in par- 
allel to an elastic spring (/io)- This viscoelastic element describes a solid with a single relaxation 
time. 




Fig. 2 Sketch of the generalized Maxwell model. The spring /io represents the shear modulus at 
high temperature. 



used the strain relations Ers = Ors/i^rs ^nd Enl = (^nl/V- After performing the 
Fourier transformation we obtain 



i(oe = -1(00/11^^ + 0/7]. 



(2) 



It follows that the shear modulus of the combined system of such a Debye relaxor 
(DR) is jUdr = £/o. 
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with relaxation time Tdr = 77 /^rs- When the viscoelastic material exhibits a sin- 
gle dominating relaxation time, then it is sufficient to consider the whole solid 



as a parallel connection between the elastic part connected in parallel to a De- 
bye relaxor, see Fig.[Tlb). Since the shear modulus is additive when connected in 
parallel, the total shear modulus becomes 



li{0)) =^o + Mdr(«)=Ao 
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1 - /COTdr 



(4) 



with g = /iRs/Ao and /Jq = /io + /^rs is the dressed elastic shear modulus. In the 
torsion oscillator and shear modulus experiments the dissipation peak is usually 
broader than that obtained from a single Debye relaxor. Thus we use a distribution 
of relaxation times attributed to a distribution of viscoelastic components with 
their own properties. To consider the general case, we need to consider a series 
of Debye relaxors with different relaxation times connected in parallel as shown 
in Fig. 12] Therefore the total contribution from the anelastic part of n constituents 

in series is given by jj.ae = LnM^"' = MrsL«[1 — 1/(1 ^ '^'^dr)]- ^he continuous 
version of this expression, when assuming a distribution of relaxation times, P{t), 
is then 



/•oo 

jU„^(a)) = jURs / dtP{ 
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(5) 



To make progress, we consider a specific form for P{t). We take the Cole-Cole 
distribution for dielectric glassesi^: 



P(t) 
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l+f2a + 2f«cosa7r' 
In this case, the contribution from the anelastic part to the shear modulus is simply 

Mrs 



^iae{0)) = - 



(l(OT) 



(7) 



Finally, the total shear modulus of the system is essentially that of the anelastic 
component connected in parallel to the elastic component of the solid. 



^i{co) = jQo 
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1- (iCOTY 



(8) 



Expression ^ is identical to the one obtained from a linear-response theory in 
which an overdamped glass component is treated as a back-action term.— This is 
not a coincidence, since the back-action term describes the response of glass to 
the external stress and glass is a anelastic material. 



3 Relaxation dynamics 

The expression obtained for the shear modulus in the previous section describes 
the dynamics of the system by one single scaling parameter x = COT. Hence all 
response quantities described by x are universal and independent of applied fre- 

22 

quency . 



We now discuss response properties that follow directly from Eq. ([Hi. The ex- 
perimental observables are the amplitude of the shear modulus, |/i |, and the phase 
delay between the input and read-out signal, (/) = arg (/j); </> measures the dissipa- 
tion of the system, which is related to the inverse of the quality factor Q~^ = tan (/) . 
Defining Ajj. as the change in shear modulus amplitude between x = and x — !> oo 
one has 

^ = \^i{x^oo)-^{x = 0)\/fio=g, (9) 

Mo 

where g measures the strength of the back-action as well as the number of con- 
stituents. At finite frequency (0 > 0, x ^ °o corresponds to the low temperature 
limit, whereas x — > is the high temperature limit. Thus the applied frequency 
defines the meaning of "high" and "low" temperature regime. 

Finally, we can estimate the height of the dissipation peak. The dissipation 
peak is centered at x = 1 and vanishes at x = and x — )• oo. Hence the peak height 



A^ = ^{x=l)-^{x = 0)=argn-j^\. (10) 

For g -^l, after some straightforward algebra, it reduces to 

40 = -^ cot(a7r/4). (11) 

For 1 < a < 2, this simplifies the expression further to 

4(/>«^cot(a7r/4)«(l-a/2)('^V (12) 

where z\0 is in units of radians. The peak height A^ depends only on the model 
parameters a and g. Note that both Ajj. and A^ are independent of frequency. 

Another important outcome of having a single scaling parameter x in the prob- 
lem is that the universal Cole-Cole plot is independent of applied frequency. We 
will elaborate on this important behavior and the associated consequences in the 
next section. 



4 Results 

We compare our theoretical calculations with the experimental data obtained by 
Day et al^o at applied frequencies of 2000 Hz, 200 Hz, 20 Hz, 2 Hz and 0.5 Hz. In 
particular we like to show the calculation utilizing Vogel-Fulcher-Tamman (VFT) 
type of relaxation: 

Here To is the attempt time, A is the activation energy, and 7b is the bare glass 
transition temperature. The VFT form represents a thermally activated process of 
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Fig. 3 (Color online) Experimental data and theoretical calculations of the shear modulus 
vs. temperature assuming a VFT relaxation time. The red and blue squares are the experi- 
mental data for the amplitude and dissipation of shear modulus. The black-solid lines show 
the theoretical calculation. The theoretical calculations use the set of parameters a = 1.31, 
g = 1.44 X 10-', /io = 0.47 pA Hz"', Tq = 50.0 ns, A = 1.92K, and Tq = -69.3 mK. Notice a 
negative To means that there is no true phase transition occurring at finite temperatures, probably 
because of strong quantum fluctuations of helium atoms. The shear modulus amplitude and the 
phase angle plotted are shifted by 0. 1 pA Hz^ ' and 3 degree with respect to the 2 Hz data. 



relaxors. The results are shown in Fig.|3] The data for 2000 Hz and 0.5 Hz show 
qualitative differences from the other three consistent data sets (200 Hz, 20 Hz 
and 2 Hz)^. We use the parameter set that best describe the 200 Hz, 20 Hz and 
2 Hz data and plot the outcome from the same parameter set for 2000 Hz and 0.5 
Hz just for comparison. 

We define the crossover temperature Tj- as the temperature where the dissi- 
pation peaks. The calculated Tx is slightly higher than in experiment for the 200 
Hz datasets, while slightly below the 20 Hz and 2 Hz dataset. As expected Tx 
decreases with decreasing applied frequency. 

Notice that although we show the results for VFT relaxation, it is not the only 
possible relaxation process that can describe the current published data; power-law 
or other types of relaxation can give similar level of agreement to the experiment 
(taking into account the number of parameters and the number of data sets be- 
ing described). In fact, a phenomenological relaxation formulation based also on 
an anelastic model was used by Syshchenko et al.^". With a significantly broad 
logarithmic-normal distribution of activation energies in the relaxation formula- 
tion, the main feature is also captured. On the other side, Iwasa proposed a relax- 
ation process ^^ based entirely on dislocation theory of Granato and Lucke^^. As 
we have discussed intensively in our previous publication '^, further experiments 



0.5 
0.4 

^ 0.3 

< 

^ 0.2 

i; 0.1 

■i 0.0 
-0.1 
-0.2 



Syshchenko et al. 2010 




• 2 Hz 20 Hz 200 Hz 




VFT calculation 




2 Hz -A- 20 Hz 200 Hz 










., 


%. 


- 


•«5 • 


■ 


• 



0.0 0.2 0.4 0.6 0.8 
Re[|//-;/J/A//] 



1.0 1.2 



Fig. 4 (Color online) The Cole-Cole plots for experimental data and for VFT calculation. For 
given form of t, all different frequency curves collapse onto one single master curve reflecting 
that (01 is the only scaling parameter. The Cole-Cole plots show reflection symmetry about 
Re[|/x — flo\/Afi]=0.5, which is a consequence of the Cole-Cole distribution function. 



at lower frequencies and lower temperature are required to characterize the exact 
type of relaxation. 

Finally, we show that the Cole-Cole plots for all three frequencies collapse 
onto one single curve, see Fig. |4l This is a consequence of x = (OZ being the 
single scaling parameter of the problem, as mentioned in Sec. 3. Note that this 
feature is universal and should be seen in other dynamic measurements besides 
shear modulus and torsional oscillator experiments.— 



5 Conclusions 

In summary, we have shown that the back-action term of a glass contribution can 
be described by a viscoelastic model through the incorporation of anelastic ele- 
ments in the constitutive equations. When considering the overdamped viscous 
property of a glass coupling to the elastic helium solid, a generalized Maxwell 
model leads to an identical expression for the complex shear modulus as the glassy 
back-action term in the linear response theory of elastic deformations. We have 
also shown that taking (OT as the scaling parameter of the problem results in a uni- 
versal scaling behavior of the real and imaginary part of the shear modulus, which 
can be seen in the Cole-Cole plot typical for glasses. Finally, we have shown that 
our calculations assuming a Vogel-Fulcher-Tammann relaxation time are in accor- 
dance with available experimental data. 
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